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The remarkable properties of silicon have made 
it the central material for the fabrication of cur- 
rent microelectronic devices. Silicon's fundamen- 
tal properties also make it an attractive option 
for the development of devices for spintronics [lj 
and quantum information processing [2145]. The 
ability to manipulate and measure spins of single 
electrons is crucial for these applications. Here 
we report the manipulation and measurement of 
a single spin in a quantum dot fabricated in a 
silicon/silicon-germanium heterostructure. We 
demonstrate that the rate of loading of electrons 
into the device can be tuned over an order of mag- 
nitude using a gate voltage, that the spin state 
of the loaded electron depends systematically on 
the loading voltage level, and that this tunabil- 
ity arises because electron spins can be loaded 
through excited orbital states of the quantum dot. 
The longitudinal spin relaxation time T\ is mea- 
sured using single-shot pulsed techniques [6j and 
found to be ~ 3 seconds at a field of 1.85 Tesla. 
The demonstration of single spin measurement as 
well as a long spin relaxation time and tunabil- 
ity of the loading are all favorable properties for 
spintronics and quantum information processing 
applications. 

Silicon is a material in which spin qubits are expected 
to have long coherence times, thanks to the predomi- 
nance of a spin-zero nuclear isotope and relatively weak 
spin-orbit coupling. However, silicon quantum dots have 
yet to demonstrate the reproducibility and controllabil- 
ity achieved in gallium arsenide devices [THE]. Here, 
we demonstrate the control and manipulation of spin 
states of single electrons in a silicon/silicon- germanium 
(Si/SiGe) quantum dot and report the first single-shot 
measurements of the longitudinal spin relaxation time 
Ti in such devices. We also show that the presence of a 
relatively low-lying spin-split orbital excited state in the 
dot can be exploited to increase the speed and tunability 
of the loading of spins into the dot. Our results demon- 
strate that Si/SiGe quantum dots can be fabricated that 
are sufficiently tunable to enable single-electron manip- 
ulation and measurement, and that long spin relaxation 
times are consistent with the orbital and/or valley exci- 
tation energies in these systems. 

The measurements we report were performed on a 
gate-defined quantum dot with the gate configuration 
shown in Fig. [T^i, tuned to be in the single-dot regime. 
The dot is measured at low temperature and in a paral- 
lel magnetic field. As shown in Fig. [TJd, an electron can 
be loaded into one of four energy eigenstates; we denote 
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FIG. 1: Device and spin-state spectroscopy, a, Scan- 
ning electron micrograph of a device identical to that used 
here. The single quantum dot location is indicated by the 
white circle. The dot barriers are tuned so that electrons 
can tunnel to and from the left reservoir only, and charge 
sensing is performed by measuring the current labeled /qpc- 
b, Voltage pulses of amplitude V pp , applied to either gate L 
or PL, adjust the energy levels of the dot. The case shown 
corresponds to the possibility of loading an electron during 
the positive phase of the pulse (left) through any of the four 
states | t d)i I i 9)1 I t e ) 5 I i e ) that are below the Fermi en- 
ergy Ef, and then unloading the electron during the negative 
phase of the pulse (right), a process that will always occur 
from I f g) or | 4- g), because of the fast orbital decay rate, c 
&; d, Measurements of gL/qpc/^Vl, in the presence of a pulsed 
voltage Vpp on gate L, as a function of VL and the amplitude 
Vpp (shown in c) or the magnetic field B (shown in d). The 
three lowest eigenstates, | \. g), | t #) 5 and | 4 e), are clearly 
visible in c, where the magnetic field B = 1.5 T. The state 
I t e ) is expected at the location shown by the red dashed 
line but is invisible due to the strong coupling of the excited 
state |e) (see Supplemental information). As a function of B 
(and with Vp P = 0.125 V), the ground and first orbital excited 
states split linearly due to the Zeeman effect. 



the states, in order of increasing energy, as | \, | t #)> 
I I e), and | t e), where the first index refers to spin 
(I having lower energy than t ) and the second to the 
ground (g) and excited (e) orbital levels. We obtain an 
experimental map of these states by measuring the differ- 
ential current dlc^pc / dVi, through a charge sensing quan- 
tum point contact while applying square voltage pulses 
to gate L. The grayscale plots of dlq^c / f dVh in Fig.[l]3,d 
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reveal dark lines when energy levels on the dot come into 
resonance with the Fermi energy [H] . The measured 
calibration factor for gate L, a = 0.125 ± 0.006 eV/V 
(see Supplemental information), translates the positions 
of these lines into a spectroscopy of the dot energy levels. 

In Fig. [TJ^, the darkest line is the orbital excited state 
\e) of energy 311 ± 19 fieV. This relatively large energy 
splitting is favorable for applications in which spin coher- 
ence is desirable [12], and it is notable because the low- 
lying excited orbital states in Si/SiGe quantum dots have 
fundamental differences compared to those in GaAs [13] , 
because of the role of valley degrees of freedom [T4H20] . 
The Zeeman splitting of the ground and excited states is 
shown in Fig. [T]i. 

A main focus of this work is to demonstrate that orbital 
excited states can be exploited to control both the rate of 
loading and the spin state of the electron loaded into the 
dot. To measure the spin state of individual electrons, 
we use a pulsed-gate technique, pioneered in Ref. [6] and 
summarized in Fig. [2^i. An electron is loaded into either 
the spin-| or spin-| ground states, | I g) or | g). Spin 
readout is performed by shifting only | t g) above the 
Fermi level, so that a spin-^ electron will tunnel out of 
the dot, but a spin-| electron will remain. When a spin-t 
exits the dot, a spin-^ will reload, because | 4- g) is below 
Ey- The readout signal is a pulse in 2qpc that persists 
as long as the dot is unloaded. 

We show in Fig. [2}d, c that loading also can be per- 
formed into the spin-split excited states, | | e) and | t e), 
and that loading into these states is much faster than 
loading into the ground state. The loading rate can be 
determined by averaging many loading events together, 
resulting in an exponential decay of /qpc from a magni- 
tude corresponding to an unloaded electron to that cor- 
responding to a loaded electron. Panels d-f show that 
the loading rates into the excited orbital states | | e) and 
| f e) are much faster than that into the spin-| ground 
state | t g)- 

The rate at which an electron loads into the dot, as well 
as the spin state of that electron, both depend strongly 
on the voltage at which the electron is loaded. In Fig. [3] 
we show that these features are observed over a wide 
range in gate voltage. Fig. [3^i shows that the loading 
rate ri oa d, measured in the same way as in Fig.[2]i-f, but 
now as a continuous function of loading level, increases 
by over an order of magnitude when the electron is loaded 
through an orbital excited state instead of through the 
orbital ground state. The loading rate shows two strong 
peaks as a function of loading voltage, corresponding to 
the excited states | t e) and | | e). Because orbital 
relaxation is very fast compared to the experimental time 
scale [2TJ [22] , one expects any electron loaded into a spin- 
t state to be measured in the state | j" g} : and any electron 
loaded into a spin-| state is measured in the state \ I g). 

The relative fraction of spin-| and spin-| loaded can be 
tuned, because the loading rate for each spin type is volt- 
age dependent. We demonstrate this directly by using 
single-shot spin readout to measure the spin-^ fraction 
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FIG. 2: Single-shot spin readout and tunable loading 
rates. Measurement of the spin state of an electron loaded 
into the dot is performed by applying a 3-level pulse sequence 
to gate PL [6]. Starting with the electron unloaded, Vpl is 
rapidly changed to a load level. Panels a, b, c correspond 
respectively to preferential loading of states | f g) , | \. e) , and 
| t e) at load voltages Vi oa d = +175, +325, and +425 mV. 
Vpl is then changed to the readout level (+75 mV), where the 
energy of a spin-t electron is above E-p, whereas the energy 
of a spin-| electron is below Ef. During the readout stage, a 
spin-t electron will tunnel off the dot, and a spin-^ electron 
will tunnel on to replace it, resulting in the current pulse 
shown in the readout phase for a and c. In contrast, a spin-^ 
electron will remain in the dot and no current pulse will occur 
during the readout phase, as shown in b. During the empty 
stage, the electron tunnels off the dot regardless of its spin 
orientation. For all cases: ti oa d = 30 ms, £ rea dout = 200 ms, 
tempty = 50 ms, and V em pt y — —200 mV. d, e, f Average of 
500 time traces of electron loading events for the load voltages 
of panels a, b, and c, respectively. Loading at the excited 
orbital state (shown in both e and f, with tunnel rates of 
1186 ± 3 Hz and 958 ± 2 Hz, respectively), is much faster than 
loading at the ground state (shown in d, with a tunnel rate 
of 106.7 ± 0.2 Hz). The tunnel rates and uncertainties come 
from exponential fits. 



as a function of the loading level, as shown in Fig. [3}d. 
The fraction of spin-^ electrons has two clear peaks as 
a function of the loading voltage, one corresponding to 
| t 9) an d a second to | t e). These peaks arise be- 
cause loading rates into specific spin states in the dot are 
a maximum when the state is near resonance with the 
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Fermi level. 

The variation in both the total loading rate and the 
spin-t fraction can be understood by calculating the load- 
ing rate for each spin state. The probability of loading a 
spin-t electron is r^/Ti ad 5 where T^ is the total rate for 
all spin-t channels, and Ti oa d is the sum of the rates for 
all ways to load the dot. A global fit to the data in terms 
of loading through spin-split ground and excited states 
with Lorentzian broadening is shown in Fig. [3^i and b. 
The fit is to both the total loading rate and the fraction 
of spin-t electrons. We treat the loading rate contribu- 
tion from each state as a convolution of a Fermi-Dirac 
distribution, a Lorentzian lineshape, and a linearized en- 
ergy dependent tunneling function [23 : 
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where E{ is the linearized energy dependent tunneling 
coefficient that relates to the transparency of the bar- 
rier [24]; 7^ is the full width at half maximum of the 
Lorentzian; Tq^ is the amplitude; V is the loading volt- 
age; and SVi = V — V^, where V{ is the position of the 
state (see Supplemental information for additional de- 
tails). The fit is in good agreement with the experimen- 
tal measurements for the loading rates and the fractions 
of spin-t and spin-| electrons, supporting the interpreta- 
tion of the data in terms of loading through specific spin 
states of the dot. 

The black line in Fig. (3}d is corrected for the finite 
spin lifetime T\ of the electrons loaded into the dot. We 
determine T\ by measuring the spin-^ fraction as a func- 
tion of the duration of the loading pulse, t\ oa ^ [6]. The 
inset of Fig. ^jp shows a typical result for a magnetic 
field of 1.85 T, with an exponential fit yielding the value 
Ti = 2.8 ± 0.3 s. For these data, V[ oad = +425 mV, so 
that loading occurs predominantly through the excited 
orbital states. We have measured T\ at different load- 
ing voltages and also measured T\ for a four-step pulse 
sequence in which the dot is held at a fourth "wait" volt- 
age. The values of T\ were all of order seconds, except 
when both the "load" and "wait" voltages were set so the 
state | t g) was aligned just below Ep, where the value 
of T\ is reproducibly much shorter: T\ = 136 ± 22 ms. 
In principle, virtual hopping of electrons from the dot 
to the leads can limit Ti [25 , but the predicted magni- 
tude is much too small to explain the observed effect, and 
this mechanism is not consistent with the behavior when 
the "wait" voltage is varied. The possibility that the 
excited orbital states are long-lived is ruled out by per- 
forming single-shot measurement with the readout level 
positioned such that | | g) and | t 9) are below Ep and 
| \, e) and | t e ) are above Ep. In this situation, even 
for short loading times (ti oa d = 20 ms), less than 0.5% of 
traces show events. 

It has been observed that T\ can be limited by dipolar 
coupling to nearby spins [26] ; we believe that the shorter 
T\ time measured when the dot is loaded through the 
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FIG. 3: Spin selective loading and spin lifetime mea- 
surement, a Plot of electron loading rate Tioad versus load- 
ing voltage Vioad on gate PL at B = 1.85 T. Error bars are 
the standard deviation of four measurements. Peaks occur 
when levels in the dot are made available for loading, en- 
abling tuning of both the loaded-spin fraction and the load- 
ing rate. The spin-split excited orbital levels | f e ) an d 
| 4- e ) load ~ 10 times faster than the ground state. The 
black solid line in a is a fit obtained by treating the loading 
rate contribution from each state as a convolution of a Fermi- 
Dirac distribution, Lorentzian lineshape, and a linearized en- 
ergy dependent tunneling function, b Spin-up fraction versus 
loading level measured using single-shot readout (see Fig. [5]), 
with tioad — 100 ms, £ rea dout — 200 ms, £ em pty — 50 ms, 
Vreadout = +85 mV and Kmpty = —200 mV. Each data point 
corresponds to an analysis of 1000 single-shot traces. The 
red dashed line is a fit using the same parameters as in a; 
the black solid line is the same fit with corrections for the T\ 
decay before measurement, a measurement fidelity less than 
unity, and a non-zero dark count rate. Error bars are y/~M, 
where M is the number of readout events, and details of the 
fitting can be found in the Supplemental information. Inset: 
Spin-up fraction versus loading time with Vioad = +425 mV. 
The solid line is an exponential fit that yields Ti = 2.8 + 0.3s. 



state | t d) is due to interaction with a nearby spin trap, 
which is suppressed by the deeper pulse associated with 
loading at the excited orbital states. The ability to load 
an electron into the dot through the excited state thus 
provides an immediate benefit, in that the additional tun- 
ability enables one to avoid sample imperfections that 
lead to a shorter spin relaxation time. 

We have demonstrated the ability to manipulate and 
measure the spin state of individual electrons in a Si/SiGe 
quantum dot, and report the first single-shot measure- 
ments of the longitudinal spin relaxation time Ti in such 
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devices. The stability and tunability of the dot enable 
tunable loading rates and tunable spin loading. These 
capabilities enhance the prospects for the development 
of Si/SiGe devices for spintronics and quantum informa- 
tion processing applications. 

Methods Summary 

Fabrication method: The dot was fabricated from 
a Si/SiGe heterostructure using the methods outlined 
in Ref. [27]. Magnetoresistance measurements using a 
Hall bar fabricated on the same heterostructure gave 
a carrier density of 4 x 10 -11 cm -2 and a mobility of 
50,000 cm 2 /Vs. The Pd top-gates were defined by elec- 
tron beam lithography in the pattern shown in Fig. [T^i. 
Measurements were performed in a dilution refrigerator 
with a base temperature of ~ 20 mK. The electron tem- 
perature was 143 ± 10 mK, measured by fitting the width 
of electronic transitions as the base temperature of the 
dilution refrigerator was raised. The dot was tuned to the 
few-electron regime, and we estimate that we measured 
transitions between two-and three-electron occupation of 
the dot. Two of the electrons formed a "closed shell" [28] 
and were essentially inert, and our measurements were 
sensitive to the state of the third electron moving into 
and out of the dot. 

Experimental method: Pulsed gate voltages were sup- 
plied using a Tektronix AFG3022B arbitrary waveform 
generator. The tunnel barriers of the dot were tuned 
so that the tunnel rates are in the range 10-1000 Hz, 
substantially slower than the inverse of the rise time of 
the pulse generator, which was < 18 ns. The time res- 



olution of the data acquisition was approximately 1 ms; 
it was limited by the bandwidth of the current preamp 
(DL Instruments Model 1211), which in turn was lim- 
ited by the need to avoid ringing in the amplification 
circuit. While the spin relaxation time can be detected 
using time-averaged measurements [29] [30], single-shot 
measurement provides a direct measurement of the spin- 
t fraction. Identification of loading and unloading events 
during the single-shot readout phase was performed by 
wavelet analysis (see Supplemental information for de- 
tails and for comparison to thresholding analysis). The 
charge sensing signal from the low-frequency quantum 
point contact (QPC) corresponds to ~ 1 x 10 -3 e 2 /h. 
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Supplemental Information 

I. ANALYSIS OF REAL-TIME ELECTRON COUNTING DATA. 

All of the single-shot counting experiments in the paper required detecting single electron tunneling events, as 
indicated by changes in the quantum point contact (QPC) current /qpc- Two different algorithms were employed 
for this determination. The first algorithm used wavelets to determine transition edges and, subsequently, levels in 
the QPC signal. All of the data in the paper were analyzed with this wavelet algorithm. In the second algorithm, 
the data were inspected using a thresholding technique. Instead of detecting edges of transitions in the signal, this 
method compares the value of the QPC signal to a threshold value and thereby infers the state of the dot. Both 
algorithms are described in greater detail below, and we show that the results of the two approaches to the analysis 
are consistent. 



A. Electron tunnel event detection using wavelets 

A single electron tunneling onto or off of the quantum dot causes changes in the magnitude of /qpc> leading to 
steps or "edges" in /qpc(^)- Wavelet analysis is a useful tool to determine the significance of edges. In fact, the 
classic Canny edge detection algorithm [31 can be interpreted in the context of wavelets. Based on observation of the 
real wavelet transform of our data and its noise characteristics, a robust edge detection algorithm was constructed 
based on the Stanford Wavelab library [32]. The algorithm has strong similarities to the Canny algorithm. First 
derivative of Gaussian wavelets [33 are used, implemented under the name 'DerGauss' in Wavelab. Significant edges 
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FIG. SI: Single-shot trace analysis using wavelets and thresholding, a &; b, Raw data and reconstructed signals 
(shown in black) using the real wavelet transform for the two single-shot traces presented in Fig. 2a, b of the main text, c, 
Real wavelet transform and modulus maxima (shown in white) for the trace shown in a. d, Results of the T\ measurement 
presented in the inset to Fig. 3b of the main text, analysed using both wavelet and thresholding algorithms. Exponential fits 
yield T\ — 2.75 =b 0.33 s for the wavelet analysis results and T± = 2.72 ± 0.37 s for the thresholding analysis results. 



are indicated by (i) the presence of a modulus maxima track over many scales of the wavelet transform, and (ii) the 
amplitude of the wavelet transform (Aw) at the location of its track [33] . 

The design of the algorithm and the tuning of its parameters were initially accomplished by the selection of a set 
of single-shot traces displaying a wide variety of inputs. These inputs included traces with events on short, middle, 
and long time scales; traces with no, few, or several events; and traces with a variety of significant features that did 
not correspond to actual events. The wavelet transforms of all traces in this set were observed, which allowed the 
construction of the algorithm and the initial tuning of detection parameters. Subsequently, a few thousand traces 
were analyzed and verified by hand, and parameters were tweaked to approach a detection rate of 98%. It should be 
noted that, in all steps of the algorithm but the last, we try to err in the direction of selecting too many (possibly 
spurious) edges, waiting until the last step to reject spurious events. 

The sampling rate employed in this experiment was 40 kHz, much higher than the 1 ms rise time of the current 
amplifier. To speed up processing, bins of ten data points were averaged, which did not result in a loss of significant 
features. The 'DerGauss' real wavelet transform was computed using four scales per octave, which provides adequate 
resolution in the wavelet domain to make real features stand out visually. The modulus maxima of this wavelet 
transformation is then computed. A sample plot of a wavelet transformation with its modulus maxima superimposed 



can be seen in Fig. SI z for the trace shown in a. 

The algorithm then iterates from the coarsest to the finest scale of the wavelet transform. At each scale, the wavelet 
amplitude squared (Ayy) for each modulus maxima track is computed. Tracks for which this value reaches ten times 
the median value of all tracks at that scale are stored for further processing. The purpose of this step is to select edges 
that are significantly larger than all other edges at the time scale under examination, which is reasonable because 
few events are expected in any particular trace. For the results presented in this paper, increasing or decreasing the 
cutoff value by a factor of two results in less than a 5% change in the correct detection ratio. 

The stored tracks are then followed to the next finer time scale, a procedure which is repeated until the best estimate 
of position for an edge is obtained at the finest time scale. For each modulus maxima track this is accomplished by 
finding the closest track for which the wavelet amplitude has the same sign on both scales. Some tracks representing 
spurious features, such as an insignificant change in the base level of the signal, or the remainder of a crosstalk edge at 
the beginning or end of the data, may fall off the edge of the wavelet transform. For this reason, all tracks that come 
to within 2.5 rise times of the edge are rejected. If cross talk from the applied pulse sequence were not significant, 
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this parameter could be set much less aggressively, as in that case only tracks that truly fall off the edge need to be 
rejected. Halving or doubling this parameter causes up to a 10% change in correct detections. At each time scale in 
the iteration, new tracks that suddenly reach the cutoff discussed in the previous paragraph are stored and considered 
in all subsequent iterations. For each track followed, is summed from the scale at which it is detected until the 
finest scale. 

When the finest time scale is reached, any edges for which the cumulative does not meet a secondary cutoff 
are rejected. This cutoff is set by finding the typical maximum value of A^ for a feature of interest by hand for a few 
traces, and requiring a sum of at least four times this amplitude. This corresponds to a significant presence of the 
edge over at least one octave. This will accept, for example, short tracks that are very strong, such as sharp spikes in 
the signal, but also longer tracks that are a bit weaker, such as long pulses with more gently sloping edges. Halving 
or doubling this parameter causes approximately a 10% change in the correct detection ratio 

All cutoffs up to this point were set to risk acceptance of edges of spurious features rather than to risk rejection of 
edges of weaker, true features. The last step in the algorithm calculates the levels between edges found, and uses this 
information to perform a final rejection based on the physics of interest. This step can have a large impact on which 
features are detected. While the algorithm used for this last step works well for the data in this paper, a completely 
different approach may be necessary for data with vastly different noise or signal characteristics. For these particular 
experiments, edges that do not border significant changes in 2qpc are rejected. This is accomplished by rejecting 
level changes that are less than five times the standard deviation of the signal between all edges. This standard 
deviation is a good measure of high frequency noise present in the signal. The sensitivity to this parameter is rather 
large, resulting in up to a 30% change in correct detection ratio when the parameter is increased or decreased by one 
standard deviation, but once set, it is robust to the variation in high frequency noise encountered in the extensive 
quantity of data analyzed for the current experiments. 

Figure [STfc and b show results from using this wavelet edge detection algorithm for the readout phase of the two 
single-shot traces presented in Fig. 2a, b of the main text. The reconstructed signals using the real wavelet transform 
are shown in black superimposed on the data. An up-edge and a down-edge are detected in a, corresponding to the 
unloading of a spin-^ electron followed by the reloading of a spin-^ electron, and no edges are detected in b. 



B. Electron tunnel event detection using thresholding 



In this algorithm, the QPC signal is compared to a threshold. A tunneling event is detected if /qpc crosses the 
threshold. If the signal remains below the threshold, zero tunneling events are assumed. While simple in concept, the 
thresholding analysis is complicated by the need to account for signal drift and high frequency noise. 

First, the data are corrected for slow drifting of the QPC signal by subtracting the value of 2qpc corresponding to 
the electron located on the quantum dot (loaded) from every trace in a given experimental run. The value to subtract 
from each trace is determined in one of two ways, depending on the variance of the signal. If the variance is small (as 
for traces with zero tunnel events), then the loaded QPC signal is determined as the mean of the data. If the variance 
is large (as for traces containing tunneling events), then a histogram of the data is produced. This histogram contains 
a noise-broadened double peak, with the peaks centered on the values of /qpc corresponding to the electron on or off 
the quantum dot. Least-squares fitting of a sum of two Guassians identifies the centers of the peaks, and the lower of 
the two center values is taken as the loaded 2qpc value. 

After the drift correction, the threshold value can be determined. A combined histogram of the QPC signal for 
many single-shot traces is produced. This histogram contains a double peak where, again, the two peaks are centered 
on the loaded and unloaded values of /qpc- The histogram is fit to a sum of two Guassian peaks, and the threshold 
value is chosen as the location of the minimum between the fitted peaks. 

Once the threshold is established, every single-shot trace can be identified as containing tunnel events (the signal 
crosses the threshold), or containing no tunnel events (the signal remains below the threshold). To reduce the 
probability that short, noise-induced spikes in the signal are misidentified as real tunnel events, only above-threshold 
signals that last longer than twice the time constant of measurement amplifier are counted. 



C. Comparison of wavelet and thresholding analysis 



All relevant data were analyzed employing both wavelet and thresholding methods. The results show good quantita- 
tive agreement (see Figure SI i). The advantage of the wavelet algorithm is that it allows each trace to be interpreted 
independently of all other traces. It is therefore uniquely suited for live, real-time analysis. 



7 




FIG. S2: Measurement of a and electron temperature, a, A charge transition peak measured with fridge temperatures 
ranging from 24 mK to 335 mK. b, The width of the charge transitions for each of the temperature traces shown in a determined 
from a global fit. Error bars for each data point are determined from the uncertainties in the fit. 



II. MEASUREMENT OF THE GATE VOLTAGE LEVER ARMS, THE ELECTRON TEMPERATURE 

AND THE ZEEMAN ENERGY 



The spectroscopy in Fig. lc,d of the main text show the electron response to an 800 kHz square wave pulse applied 
to gate L (see [IT] for details of the spectroscopy technique). The rate for an electron to tunnel onto the quantum dot 
is tuned to be slow compared to this pulse frequency. A lock-in signal with a carrier frequency of 97 Hz applied to 
gate L is used to measure diQpc/dVL in response to the square wave pulse. The dark lines indicate when additional 
energy levels become accessible resulting in a sudden increase in the electron tunneling rate. In both c and d, the 
dark contrast of the excited orbital state indicates that it is more strongly coupled to the reservoir than the ground 
orbital state, because it causes a larger increase in the electron tunnel rate and thus the electron response. Only the 
lower energy of the strongly coupled orbital excited states, | ^ e), is clearly visible in this spectroscopy, because once 
this state is accessible the tunnel rate is comparable to the pulse frequency and additional energy levels have no effect 
on the electron response. 

In order to translate the spectroscopy data into energy levels on the dot, it is necessary to determine the lever 
arm, a, for gate L. The width of a charge transition peak in gate voltage is determined by both a and the electron 
temperature [34]. In order to determine these variables independently, we fit the width of a charge transition as a 
function of temperature by increasing the temperature of the dilution refrigerator. The results for fridge temperatures 



(^fridge) ranging from 24 mK (base temperature) to 335 mK are plotted in Fig. S2i. The peak indicating the charge 



transition clearly decreases in amplitude and increases in width as the temperature is increased. 



We extract the width of each peak in Fig. S2i in terms of T e /a by performing a global fit to the data. Each peak 
is fit with the equation, 

B „„^-2 fa(V -V L ] 



dI QPC /dV h (V L ) = ^coah- [-^j^ ) + D (SI) 

where B, D and Vb are fit parameters that are the same for all the peaks, and T e /a is a fit parameter that is optimized 
for each peak independently. The fit results for T e /a for each of the peaks in Fig. [S2^ i are shown in b. As expected, 
at high temperatures T e /a is linear in Tf r id ge ? while at low temperatures T e /a saturates. We fit this cross-over using 

S2p, and the results for the 



the functional form T e = y Tq + T f ^ idge . This fit is shown as the solid blue line in Fig. 

two fit parameters, a and To give a = 0.125 ± 0.006 eV/V for gate L and T e = 143 ± 10 mK at base temperature. 

With a measure of the electron temperature, it is possible to determine the lever arm for additional gates by 
performing a single fit. We used gate PL to sweep through the same charge transition at base temperature, and the 
resulting fit gives 0.00144 ± 0.00013 eV / V for a multiplied by the attenuation (note: the voltage on gate PL is 
attenuated by a factor of ~ 30). 

Using the above results for a for gate L, it is possible to extract the Zeeman energy as a function of magnetic field, 
£?, for the ground orbital state from the excited state spectroscopy data presented in Fig. Id of the main text. We 
take vertical line cuts to determine the separation, AVl, between the peaks corresponding to the states | j g) and 



| t g) for B w 0.7 to 2 T. The results, after correcting for the gate voltage lever arm, are plotted in Fig. S3 

Similar to the results for a GaAs quantum dot presented in [35] , at high magnetic fields we observe a clear deviation 
from the bulk silicon g-factor of 2 (black dashed line) . Following the method in [35] , we fit the Zeeman splitting data 
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B(T) 



FIG. S3: Zeeman splitting versus B-field. Zeeman splitting as a function of magnetic field determined from the data 
presented in Fig. Id of the main text. Error bars are due to a combination of uncertainty in the peak separation, AVl, and 
uncertainty in the lever arm, a for gate L. The dashed black line indicates g — 2 and the solid blue line is the result of a fit 
using a second-order polynomial. 

with a second-order polynomial restricted to pass through the origin. The result gives 

g = (1.96 ± 0.15) - (0.25 ± 0.05)£(T). (S2) 
The -independent contribution of this result is in good agreement with the value for an electron spin in bulk silicon. 

III. GLOBAL FIT OF LOADING RATES AND SPIN-t FRACTION 

In this section, we give details of the modeling of the loading rates and the spin-^ fraction as a function of the 
loading gate voltage. 



A. Modeling of the loading rates 

We perform a global fit to the data for loading rate and the spin-f fraction as functions of the loading level, shown in 
Fig. 3a and b of the main paper. Fig. 3a shows two pairs of peaks that correspond to the loading of the Zeeman-split 
ground | t #)> I i 9) an d excited | \ e), \ I e) orbital states, as well as a shallow minimum on the right that we believe 
is due to a contribution from the loading of a higher orbital | t h), | I h). Fig. 3b shows two peaks as a function 



of loading level that correspond to the measurement of an increased fraction of spin-^ electrons. Tables SI and S2 
provide the numerical results and parameters for the global fit. 

In fitting to the loading rate data of Fig. 3a, we take into account energy-dependent tunneling across the barrier [24] 
as well as the lifetime broadening of energy states of the dot, and we use the Landauer-Buttiker formalism [23] to 
model the loading rate. The lifetime broadening of energy levels is modeled using a Lorentzian function. The overall 
loading rate ri oa d as a function of loading level V is given by T\ Q ^(V) = IMV), where the sum is over all possible 
transport channels i G {|, t; 9, e ? h}. 

For each of these channels, the loading rate is given by the convolution of the Fermi-Dirac distribution, the energy 
dependent tunneling function, and the Lorentzian function. It is given by 

where Ei is the linearized energy dependent tunneling coefficient that relates to the transparency of the barrier [24] ; 7^ 
is the full width at half maximum of the Lorentzian; Toi is the amplitude; V is the loading voltage; and SVi = V — V*, 
where Vi is the position of the state. 

In the fitting, the Zeeman splitting for each pair of states was identical, and initial guesses were made for the overall 
scale, position and Lorentzian width for the higher orbital states (| f h), \ ih)), with the leveling off of ri oa d at about 
600 mV as a guide. These contributions are plotted as dotted lines in Fig. 3a of the main text. After each fitting run, 
the parameter estimates for | ^ h) and | ^ h) were revised such that the ratio of overall scaling of the rates and the fit 
repeated so that self-consistency is achieved. 
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Following the observations in Ref. [10 j, the amplitudes of the loading rates for different spin states were allowed to 
differ; we find, for each orbital, that I\ < T^. 

We note that, while the energy dependent coefficients for the ground orbital states gave a good fit when they were 
constrained to be equal, this was not the case for the excited orbital states. It was necessary to have E± e > E^ e , as 
can be seen from the fact that the downward slope on the right of the highest peaks, which is strongly dependent on 
E^ e , is much steeper than those of the ground orbital states. 

The ground orbital states were modeled as delta functions, while the excited orbital states were modeled with a 
finite energy width as fit parameters, as consistent with the hypothesis of short-lived excited states. However, due to 
the fact that the temperature broadening of the Fermi function limits the resolution with which we can extract any 
energy information, only an upper bound on the energy width due to lifetime broadening could be obtained, giving a 
lifetime of at least 4 x 10 2 ps. 



B. Calculation of spin-t fraction 

In the limit of infinite loading and unloading rates, and no spin-t decay, the fraction, n^, of spin-t electrons measured 
as a function of loading level V is simply 

n t (V) = r r '° adt( ^ /v (S4) 

1 loadjtotaU v ) 

where the loading rates Ti(V) are functions of the voltage V at which the dot is loaded, r^adt is the combined rate 
for loading all spin-t channels, and ri oa d, total is the total loading rate for all spin channels. This quantity is plotted 
as the dashed red line in Fig. 3b of the main paper. 

So that the modeled spin-t fraction correctly takes into account finite rates of loading, unloading and spin-t decay, 
we set up rate equations for the load and readout phases and solve for the fraction of spin-t at the end of each phase. 

Given a pair of spin-t and spin-t states, the rate equations can be written as 



dn^(t) 
dt 

dt 



= -T x 1 n t + rioadt (1 ~ n t _ n i) - r un i oadt n t , (S5) 

= Tf^t + r^adt (1 - % - n±) ~ r un load|^. (S6) 



where n^(t) and n±(t) are the fractions of spin-t and spin-t respectively. They are explicit functions of time, and 
implicit functions of loading voltage V via the loading rates Ti(V). 

With the dot empty initially, we impose initial conditions n^(0) = ^(0) = and solve to obtain the solution, which 
is a complicated function of the rates and time. However, during the load phase, both the spin-t and spin-t states are 
below the Fermi level of the lead, so the unloading rates are much smaller than the loading rates and can be neglected: 
i.e., r un ioadt = r u nioad| = 0- This assumption leads to a simplified solution that offers more physical insight, given by 

n t (t) = _ e -tTf 1 (1 _ e -(rw t +r,o^-T 1 - 1 )t )> (S7) 

rioadt + rioadt — ^1 

(rwt - rr 1 ) (i - e -*dwiw) + r (x _ e -tr-i\ 

= = -= ==i ^ '-■ (S8) 

1 loadt + 1 loadi — 1 1 

It is clear that, in the limit of loading rates fast compared to the decay rate, the approximation w 

contributes 
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-^rdecay (r^adt/r load, total) holds. If the rates are comparable, then the last term in brackets in Eq. 
significantly. 

We can generalize to include the loading of other spin states whose population comes from the tunneling in of 
electrons from the lead and not from the decay of the spin-^ state in question. This is necessary as the spin-^ states in 
question are each sandwiched between two spin-^ states and a higher orbital state. The solution to the spin-^ fraction 
then includes the total tunneling rate into all states, r to tai = T^ad aiu + Tioad aiit? gi yen by 

n t (t) = e-^ 1 (1 - e-^-T-)*). (S9) 

J- total — J- 1 



Next, we consider the correction to the spin-^ fraction during the readout phase. During this phase, because the 
spin-t state is above the Fermi level and the spin-^ state is below it, we neglect the spin-t loading rate as well as the 
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FIG. S4: Measurement of the spin-t unloading rate. Histogram of the times at which an up-edge is detected in the 
readout phase of the single-shot traces for the data presented in Fig. 3b of the main text. The exponential fit gives the rate at 
which a spin-t electron tunnels off the dot, r~^ load ^ = 21.6 ± 1.3 ms. 

spin-^ unloading rate. Because the quantity measured in the experiment is the number of spin-^ electrons that tunnel 
out to the lead n c ^, we solve the following rate equations for n c ^(£), 



dn^(t) 
dt 

dt 
dn c ^(t) 
dt 



-T x - r un i oadt n^, (S10) 
Tf 1 ^ + r loadt (l - n t - H ), (Sll) 

(S12) 



The rate r un i oac ^ at which a spin-^ electron tunnels off the dot during the readout phase is determined by fitting 
to a histogram of the times at which up-edge is detected in the single-shot traces. The results for the data presented 
in Fig. 3b of the main text are shown in Fig. S4, where '0' corresponds to the start of the readout phase. The 



exponential fit gives r un 1 load ^ = 21.6±1.3ms. 

With the initial condition of a loaded spin-| state, n^(0) = 1, we obtain the measured spin-t fraction as 

nct(t) = T J" nl r ° adT (l - c-^+^t)) (S13) 

J-l T 1 unloadt V 7 

With an independent determination of the spin-t detection fidelity and the dark count rate (92% and 3% respectively, 
as described in the next section), the measured spin-up fraction fully corrected for finite tunnel rates during the 
load and readout phases of the 3-level pulse sequence and for spin-| decay, is given by 

« t = 0.03 + 0.92 X n t (tioad) X ™ct (Readout), (S14) 

where ti oa d and ^readout are the durations for the load and readout phases (100 ms and 200 ms respectively for the 
data in Fig. 3b). This result is plotted as the bold, black line in Fig. 3b of the main paper. 

C. Single-shot spin readout visibility 

We determine the visibility of our spin readout using the method in [6]. Since our g-factor is positive, spin- 4- is the 
lower energy spin state in our case, and we therefore analogously define a as the probability that a spin-^ is detected 
as "up" and j3 as the probability that a spin-| is detected as "down", and we use corresponding definitions for the 
two contributions to /3: f3\ and /?2- The visibility of the spin readout is then given by 1 — a — j3. 

Due to the finite electron temperature, it is possible, although rare, for a spin-^ electron to tunnel off the dot during 
the readout phase even though its energy is below Ep, causing a spin-| electron to be misidentified as spin-| . This 
'dark count' is determined by the saturation value of the exponential fit in the inset to Fig. 3b in the main text. The 
result gives a = 0.031. 

During the readout phase, if a spin-| electron tunnels off the dot and is replaced by a spin-| electron at a rate 
comparable to the bandwidth of the measurement electronics, then the resulting step in /qpc will be too small to 
detect, causing a spin-^ to be misidentified as spin-^ . As in [6], we determine the spin-^ detection fidelity by applying 
the 3-level pulse sequence in reverse, and we find it to be 91.8%. This gives /3 2 = 1 — 91.8% = 0.082. 
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TABLE SI: Results of fit 



r it Parameter 




Value [Uncertainty] 


Energy dependent coefficient 


ig), \ \ g) (meVJ 


U.14o [z/J 


Energy dependent coefficient 


^ ey ^mev ) 


u.iyo [oUj 


Energy dependent coefficient 


T e) (meVJ 


U.UD4 [zoj 


Position | | g) (mev ) 




U.ioO [loj 


Position | 4- e) (meV) 




U.40D [4zJ 


Position t e) (meV) 




(J. 585 [58] 


Amplitude ^ g) (Hz) 




116 [8] 


Amplitude | t g) (Hz) 




110 [7] 


Amplitude \. e) (Hz) 




1,916 [62] 


Amplitude | t g) (Hz) 




610 [33] 


Lorentzian width for e) (meV) 


5 xlO -6 



Quantities derived from fit parameters 

Position | I g) (meV) 0.057 [9] 

Zeeman splitting (meV) 0.129 [17] 

g-factor 1.20 [16] 

Lifetime of \e) (ps) 440 



TABLE S2: Constants and independent parameters 



temperature T (mK) 143 [10] 

magnetic field (T) 1.85 

alpha x attenuation (eV / V) 0.00144 [13] 

T hg (ms) 136 [22] 

Ti, e (ms) 2.8[3] x 10 3 

Parameters for | ^ h) 

Position (meV) 1.146 

Lorentzian width (meV) 0.433 

Amplitude (Hz) 3100 

Parameters for | t h) 

Position (meV) 1.275 

Lorentzian width (meV) 0.433 

Amplitude (Hz) 310 



Due to the finite time for an spin-^ electron to tunnel off the quantum dot during the readout phase, it is possible 
that an electron that is spin-t at the start of the readout might decay to a spin-^ before it tunnels off the dot. This 
contribution to the spin readout fidelity is given by f3\ = 1/(1 + TiT un i oa( ^). From the exponential fit in Fig. S4, we 
find r~^ load ^ = 21.6 ms. This result combined with the measured value of T\ = 2.75 s gives /3i = 0.0078. 

Finally, we determine a spin-^ detection fidelity of 1 — a ~ 0.97 and a spin-^ detection fidelity of 1 — j3 ~ 0.91, 
and this amounts to a total single-shot readout visibility of 88%. This visibility is notably improved from the 65% 
reported in [6 , and the majority of the improvement is due to the long spin lifetime compared to the relatively fast 
electron tunneling time. 
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